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ABSTRACT.  One-way  phase  observations  of  GPS  satellites  at  three  ground 
stations  were  differenced  between  stations  to  form  interferometric 
observations.  These  single-difference  observations  were  differenced  between 
satellites  to  form  double-difference  observations.  Me  analyzed  the  one-way, 
single-,  and  double-difference  observations  by  least-squares  to  determine  the 
orbits  of  the  satellites.  In  most  cases  the  formal  standard  errors  of  the 
estimated  satellite  positions  w^re  of  the  order  of  several  meters  — 
equivalent  to  a  few  parts  in  10  of  the  orbit  radius.  The  actual  errors 
of  the  estimated  orbits  are  unknown,  but  we  were  able  to  test  the  day-to-day 
precisions  by  using  these  orbits  to  analyze  observations  from  an  independent 
pair  of  stations,  to  determine  the  baseline  vector  between  them.  The  results 
of  this  analysis  confirmed  the  uncertainty  estimate  of  a  few  parts  in  10  . 
Further  confirmation  was  obtained  from  a  comparison  of  orbits  estimated  from 
disjoint,  but  interleaved,  spans  of  observations. 


INTRODUCTION 


If  the  NAVSTAR  Global  Positioning  System  (GPS)  is  to  be  useful  for 
crustal  motion  monitoring,  the  orbits  of  the  GPS  sate  I  I i tes^wi I  I  need 
to  be  known  with  uncertainties  of  the  order  of  1  part  in  10  or  less. 

This  level  of  accuracy  has  not  been  achieved.  A  major  problem  has  been  the 
instability  of  the  cesium-beam  frequency  standards  which  are  employed  in 
most  of  the  present  tracking  stations. 

To  show  that  interferometric  observations  by  stations  equipped  with 
hydrogen-maser  frequency  standards  can  yield  better  accuracy,  we  have 
installed  the  dual-band  tracking  receivers  of  the  Air  Force  Geophysics 
Laboratory  at  the  NEROC  Haystack  Observatory  (HO)  in  Massachusetts,  at  the 
U.S.  Naval  Observatory  Timing  Service  Substation  (NOTSS)  in  Florida,  and  at 
Harvard  College  Observatory's  George  R.  Agassiz  Station  (GRAS)  in  Texas. 


THE  OBSERVATIONS  AND  DATA  ANALYSIS 


Observations  made  by  these  receivers  on  four  days  in  January,  1984  and  six 
days  in  August,  1984,  have  been  analyzed  to  determine  the  orbits  of  the  GPS 
satellites  then  available.  As  the  GRAS  station  was  not  installed  until  ^ 
July,  1984,  the  January  analysis  utilized  a  single  band  (LI  only)  Macrometer'^ 
tracker  (8ock  et  al.,  1984)  operated  by  Aero  Service  in  Phoenix,  Arizona  where 
a  hydrogen  maser  was  lacking  (a  cesium-beam  standard  was  used).  The  GRAS 
station  replaced  the  Arizona  station  in  the  August  analysis.  Figure  1  shows 
the  ground  tracks  for  NAVSTARs  1,3, 4, 6, 8,  and  9.  Also  shown  are  projections 
of  the  25  degree  elevation  cutoffs  for  HO,  NOTSS  and  GRAS.  (This  figure  would 
change  only  slightly  if  Phoenix  were  shown  instead  of  GRAS.)  Highlighted 
are  the  regions  of  mutual  visibility  and  the  ground  tracks  of  the  satellites 
through  these  regions.  For  illustrative  purposes  we  have  assumed  a  25  degree 
minimum-elevation  cutoff.  This  assumption  was  intended  to  account  for 
possible  obstructions  of  the  satellites;  the  receivers  are  able  to  track  to 
fifteen  degrees  above  the  horizon. 

The  obse.'vat i  ons  consisted  of  samples  of  the  phases  of  the  reconstructed 
carriers  of  the  LI  and  L2  signals  received  from  each  satellite,  each  relative 
to  the  phase  of  a  local  reference  oscillator  at  each  site.  The  time  span  of 
the  daily  observations  of  a  given  satellite  was  between  3  and  5  hours, 
depending  on  the  visibilities  of  the  satellite  from  the  three  tracking  stations 
(Figure  1).  After  combining  the  LI  and  LZ  phase  observations  at  each  epoch  in 
order  to  remove  ionospheric  effects,  we  combined  the  data  from  different 
stations  to  make  single-difference  observations  (between  stations,  and 
nonredundant)  and  double-difference  observations  (differenced  again  between 
satellites).  (See  Bock  et  al.,  these  proceedings.) 

The  orbits  of  the  GPS  satellites  were  computed  by  direct  numerical 
integration  of  the  equations  of  motion.  For  the  January  analysis  the 
WGS72  gravity  field  (Seppelin,  1972)  was  used,  and  for  August,  the  GEM  L2 
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Macrometer  is  a  registered  trademark  of  Aero  Service  Division,  Western 
Geophysical  Company  of  America  ("Aero  Service"),  which  purchased  the  business 
and  substantially  all  the  assets  of  Macrometrics,  Inc.,  on  March  8,  1984. 


field  (Lerch  et  al.,  1982),  each  to  degree  and  order  eight.  The  effect  of 
the  direct  solar  radiation  pressure  was  accounted  for  with  an  isotropic  model 
and  with  a  free  parameter  for  each  satellite.  Me  also  included,  with  a  free 
parameter,  a  force  along  the  y-axis  of  the  spacecraft. 

The  coordinates  of  our  GPS  tracking  stations  were  determined  differently 
for  the  January  and  August  analyses.  For  the  January  analysis,  we 
obtained  MGS72  coordinates  of  NGS  monuments  which  were  close  to  each  station 
<L.  0.  Hothem,  private  communication).  These  MGS72  coordinates  had  estimated 
uncertainties  of  about  0.5  m.  Me  then  used  local  survey  ties,  performed  by  NGS, 
to  obtain  positions  of  our  antennas  relative  to  the  NGS  marks.  For  the  August 
analysis  we  defined  the  origin  of  our  reference  system  by  the  Lageos-der i ved 
coordinates  of  the  McDonald  Laser  Ranging  Station  (MLRS)  (B.  0.  Tapley  et  al., 
private  communication,  1985),  located  6  km  from  the  GRAS  receiver.  Me  then 
used  Mark  III  VLBI  determinations  of  the  baseline  vectors  to  HO  and  NOTSS  from 
GRAS  <0.  S.  Robertson,  private  communication,  1984)  and  local  area  surveys  by 
NGS  to  relate  the  GPS  and  the  VLBI  antenna  positions  at  HO  and  NOTSS.  At  GRAS, 
the  vector  between  the  GPS  antenna  and  the  nearest  NGS  mark  has  not  yet  been 
accurately  surveyed  by  NGS.  However,  when  we  were  at  GRAS  installing  the  antenna 
we  made  measurements  (with  a  steel  tape)  with  respect  to  four  NGS  survey  marks 
from  which  we  have  estimated  our  antenna's  position.  He  estimate  that  this 
crude  determination  has  an  uncertainty  of  about  20  cm. 

To  fit  the  phase  observations,  we  adjusted  simultaneously  by  least  squares 
the  six  initial  conditions  and  the  coefficients  of  the  solar  radiation  pressure 
and  y-axis  force  models  for  each  satellite,  a  clock  rate  parameter  for  each 
baseline  on  each  day,  and  an  additive  phase  bias  for  each  series  of  observations 
of  each  satellite.  The  station  coordinates  were  held  fixed  in  the  solutions. 


THE  JANUARY  ANALYSIS 


In  January  we  observed  five  GPS  satellites  (NAVSTARs  I, 3, 4, 6,  and  8). 

In  the  analysis  we  included  one-way  observables,  single-differences  and 
doub I e-di fferences. 

The  estimated  uncertainties  in  satellite  orbital  positions  are  given  in 
Table  1.  For  each  satellite  we  show  the  formal  standard  error  for  the  most 
poorly  determined  Cartesian  coordinate. 

To  test  the  accuracies  of  the  derived  orbits,  we  generated  orbital 
ephemerides  by  reintegrating  from  the  initial  conditions  estimated  in  our 
analysis;  then  we  examined  the  geodetic  results  obtained  when  we  used  our 
orbits  to  process  double-difference  data  from  observations  on  a  completely 
independent  baseline.  These  observations  were^ perf ormed  by  NGS  <C.  Goad, 
private  communication,  1984)  with  Macrometer V- 1000  Interferometric 
Surveyors  on  a  13  km  baseline  in  Southern  California.  This  baseline  was 
measured  three  times:  on  January  26,  27,  and  30,  1984.  He  processed  the 
observations  from  each  day  separately  in  order  to  see  the  day-to-day 
consistency  of  the  baseline  determinations  and,  by  implication,  of  the 
orbit  determinations. 

Me  first  determined  the  California  baselines  using  the  observations  of  all 
five  satellites.  He  noticed,  however,  that  the  residuals  for  NAVSTAR  6  with 
respect  to  the  a  priori  coordinates  of  the  baseline  were  significantly 


larger  than  for  the  other  five  satellites.  Me  attributed  this  to  a 
relatively  poor  determination  of  the  NAVSTAR  6  orbit.  Our  determination  of 
of  this  satellite's  orbit  is  suspect  since  NAVSTAR  6  passes  far  to  the  west  of 
our  tracking  network  <see  Figure  1).  We  then  repeated  our  California  baseline 
determinations  without  the  observations  of  NAVSTAR  6.  Table  2  gives  the  three 
individual  estimates  of  the  baseline  vector  which  we  obtained  using  our  orbits 
for  the  remaining  satellites.  Also  shown,  for  each  day,  are  the  root-mean- 
square  (rms)  of  the  double-difference  phase  residuals  from  the  least-squares 
adjustment  of  the  baseline  vector.  The  day-to-day  consistency  of  the  baseline 
estimates  indicate  that  our  orbit  determinations  are  precise  to  within  about 
0.3  ppm.  A  precision  of  0.3  ppm  corresponds  to  about  6  m  along-track  or 
across-track  at  the  orbital  altitude  of  20,000  km  and  is  not  inconsistent  with 
the  estimated  uncertainties  of  Table  1. 


THE  AUGUST  ANALYSIS 


In  August  we  observed  NAVSTARs  1,4, 6, 8,  and  9.  Me  had  learned  from  our 
January  analysis  that  the  one-way  phase  observations  do  not  contribute 
significantly  to  the  orbit  determination.  For  this  reason  only  the  single¬ 
differences  and  double-differences  were  included  in  the  August  analysis. 

In  the  August  analysis,  we  had  no  suitable  independent  observations  to 
test  the  accuracy  of  our  orbit  determination.  In  order  to  make  at  least  a 
consistency  check  we  divided  the  August  observations  into  2  subsets:  set 
A  consisted  of  observations  on  days  215,  217,  and  219;  set  B  with 
observations  on  days  216,  218,  and  220. 

The  estimated  uncertainties  in  satellite  orbital  positions  obtained  from 
the  two  analyses  A  and  B  are  given  in  Table  3.  Again,  the  formal  standard 
error  for  the  most  poorly  determined  Cartesian  coordinate  is  shown. 

To  test  the  accuracy  of  our  estimation,  we  reintegrated  the  orbits 
using  the  initial  conditions  estimated  from  data  sets  A  and  B,  and  then 
compared  the  orbits  by  plotting  the  radial,  along-track,  and  across-track 
differences  throughout  the  time  span  of  overlap  between  the  observation  sets. 
The  maximum  of  the  difference  in  each  component,  for  each  satellite,  is  shown 
in  Table  4.  All  of  the  maximum  along-track  differences  are  between  5  and 
10  m,  whereas  the  peak  radial  and  across-track  differences  are  typically 
about  2  m. 

If  we  assume  that  the  errors  in  the  "A"  and  "B"  orbits  are  equal,  and 
recognize  that  the  differences  between  them  are  the  sums  of  peak  (not  rms) 
errors,  then  we  conclude  that  the  rms  errors  in  the  orbits  are  no  worse 
than  3-4  m  in  any  component.  These  errors,  like  the  formal  uncertainties, 
are  about  a  factor  of  two  smaller  than  the  orbital  errors  deduced  for  the 
January  analysis. 

DISCUSSION  AND  CONCL  CIS  IONS 


The  post-fit  residuals  from  a  typical  set  (day  219)  of  single-difference 
and  double-difference  observations  from  the  August  observations  are  shown  in 
Figure  2.  The  appearance  of  systematic,  common-mode  variations  in  the 
single-difference  residuals  suggests  that  the  largest  source  of  error  in  our 


analyseG  is  the  unmodeled  fluctuations  in  the  local  oscillators.  Although 
we  can't  rule  out  a  contribution  from  the  receivers  themselves,  we  have 
independent  evidence  that  much  of  the  fluctuation  is  due  to  the  variations 
in  the  relative  phases  of  the  hydrogen  maser  frequency  standards  at  the  sites. 
These  variations  should  cancel  in  double-differences.  Other  potentially 
significant  sources  of  errors  are  variations  in  the  tropospheric  path  delay, 
any  uncorrected  cycle  slips,  and  unmodeled  accelerations  of  the  satellites. 

We  conclude  that  we  have  not  yet  reached  the  accuracy  goal  of  1-in-lO^. 

With  a  more  complete  constellation  of  satellites,  and  perhaps  additional 
tracking  stations,  it  will  be  possible  to  use  only  double-difference  observations 
to  determine  the  orbits.  Since  double-differences  are  free  of  the  effects  of 
receiving-site  frequency-standard  instability,  it  should  be  possible  to 
determine  the  orbits  more  accurately  than  we  have  been  able  to  do.  In  our 
January  and  August  experiments  we  had  too  little  double-difference  data  because 
there  was  too  little  two-site,  multiple-satellite,  mutual  visibility. 
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Figure  2.  Single-difference  and  double-difference  phase  residuals  from 
the  ilavstack-Richmond  baseline,  dav  219. 


Table  2.  Results  of  baseline  determinations.  For  each  day  the  estimated 
latitude,  longitude,  and  height  of  station  MOJAVE  1  are  given, 
assuming  latitude  35  14  54.00722,  longitude  116  47  27.95387, 
and  height  1049.805  m  for  station  MOBLAS  which  was  at  the  origin 
of  the  baseline.  The  length  of  the  baseline  was  about  13  km. 


Day 

Lati tude 

Long i tude 

Height  rms  of  res 

(mm) 

26 

35°  19'  53".  99396 

116  53  23.94062 

906.568 

5.1 

27 

.99381 

.94087 

.520* 

5.2 

30 

.99389 

.94073 

.564 

7.8 

Mean 

.99389 

.94074 

.566 

6.0 

S.D. 

0.00008 
<0.18  ppm) 

0.00013 
(0.23  ppm) 

0.003 
(0.22  ppm) 

— 

his  height  excluded  from  computations  of  mean  and  standard 
iviation  (S.D.);  we  suspect  a  blunder  in  the  field  measurement  ^ 
itenna  height  above  survey  mark.  This  same  anomaly  is  seen  when 
ihemerides  from  the  Naval  Surface  Weapons  Center  are  used  to 
termine  the  baselines. 


Table  3. 


Ana  lysis 


Estimated  uncertainties  in  satellite  orbital  positions 
<one  sigma;  worst  of  three  coordinates)  from  two  analyses 
of  observations  at  HO,  NOTSS,  and  GRAS,  August  1984. 


A 
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Days  215,217,219  Days  216,218,220 
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Table  4.  Differences  in  satellite  position  within  span  of 
overlap  (days  216-213)  between  orbits  determined 
from  A  (days  215,217,219)  and  B  (days  216,218,220) 
data  sets. 
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